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BUCKLING THE ELASTO-PLASTIC RANGE 


Vincenzo 


SYNOPSIS 


safe value the buckling load for beam the plastic range 
rigorously proved given the tangent modulus theory Engesser 
rather than the theory von Karman previously 
shown Shanley for simplified beam model. 


INTRODUCTION 


The critical value the buckling load for simple beam derived 
the elastic range under the Euler assumption that the compressive 
load remains constant during buckling. Experiments prove, instead, 
that during buckling the external compressive load can increase. the 
present paper the buckling load simple beam evaluated without 
assuming the load remain constant, and proved that, 
ance with the results obtained Shanley for simplified beam model” 
that the buckling load must obtained the tangent modulus formula 
Engesser, even the buckling stresses are above the elastic limit. 


The Differential Equation Elasto-Plastic Buckling 


Consider simple beam constant cross-section under the action 
two equal compressive axial forces and let the value the 
thrusts which buckling will start. will assumed that, the 
beam buckles, the load increases and the increase given 
stage the lateral deflection will indicated AP. The stress 
the beam before buckling constant and equal where 
the area the cross-section the beam. After buckling the stresses 
due bending the side the beam will related the 
strain the law while the stresses the tension side 
the beam will equal these formulas the tangent modu- 
lus and Young’s modulus for the beam material. The bending 
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stresses have resultant along the axis the beam, constant 
along the beam length, and produce bending moment varying along 
the beam length, magnitude AP, where the lateral deflection 
the beam. Indicating (Fig. the distance from the neutral axis 
the boundary fiber whose compression reduced the bending 
stress, and assuming that plane section remain plane, may write: 


where the maximum compressive the maximum tensile 


thea 


bho, 


means Eq. (1) the bending moment may written: 


(2) 


may seen immediately that for =1/2 and h/12, 
that: 


For AP=0, E/E, and hence: 


These are the well known reduced modulus relations von Karman. 

These results hold for <h. For h<hg the are compres- 
sive all over the section and Young’s modulus does not appear the 
formulas since over the whole section. The moment instead 
being given Eq. (2), given 


(3) 


4 
= 
= 
(1) 
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where now 


The distance can take any value between and where 


the value the von Karman relations, i.e., the value cor- 
The function (hg) graphed Fig. with its two asymptotes 
hy and A = 0. 
The radius curvature section assuming plane sections 
remain plane, related the strains the equation 


from which 


Hence, Eq. (1): 


where the moment inertia the cross-section with re- 
spect the centroidal axis parallel the neutral axis. The quantity: 


has dimensions and may considered ideal elastic modu- 
lus. Here again Eq. (5) valid for while for have 
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this case seen that the modulus coincides with the reduced 
modulus von Karman. may therefore conclude that 
varies between and varies between and and that for 

The function (hg) graphed Fig. while the function (A) 
appears Fig. 

Indicating the thrust Pg+ corresponding lateral dis- 
placement v(z), may let: 

(7) 


For (z) identically equal zero and this expression for 
becomes indeterminate. Let therefore examine the conditions under 
which the beam may start buckling. Indicating infinitesimal 
increase the load the stresses due are the type given 
Fig. the end sections the b2am, where the are uni- 
formly distributed; the other sections, instead, their distribution 
linear since they give rise bending moment vXdP. But may 
easily seen that cannot positive over the whole section all sec- 
tions the beam. Consider this purpose the limiting case occuring 
when becomes zero one the extreme fibers, and split into 
two forces and dP", the first corresponding uniform distribu- 
tion equal for all sections, and the other corresponding 
with zero value one the extreme fibers. The bending moment 
any section would then given and changing the load 
contrary the assumption that deflection starts Pg. The section 
where zero one extreme fiber is, symmetry, the middle sec- 
tion the beam (Fig. 4). buckling starts with this kind distri- 
bution, every section 


and hence, every section E;. This condition gives the lower lim- 
for the buckling load the beam. 

obtain the upper limit when every section and hence 
E'= Ey. this case the stress distribution the beam sections are 
all similar. The intersections the neutral axis (Fig. all lie ona 
straight line parallel the centroidal axis distance from the ex- 
treme fiber where compressive. this limiting and 
the value indeterminate. 

All other possible cases fall between the two extreme cases consid- 
ered above. other words, buckling may start with any value the 
ideal modulus provided 
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From these results may seen that continuous, uniform 
function hg; but from Eq. (2), also uniform, continuous function 
hg: hence continuous, uniform function and, Eq. (7), 
also and the parameter must noticed that changes sign 
with while the signs and remain unchanged: hence 

means the definition curvature, valid for small strains, 


and Eq. (5), find that 


Hence the differential equation elasto-plastic buckling becomes: 


(8) 
while the boundary conditions are v(0)= and =0. Letting: 


and indicating means double prime the derivative with respect 
Eq. (8) becomes: 


(9) 
The boundary conditions for Eq. (9) are: 


z=0 v=0 
v=0 (10) 


The non-linear, second order differential equation (9) becomes Euler’s 
equation for const. 

The above given treatment the buckling problems holds for beams 
rectangular cross-section, but the results remain substantially the 
same for beams arbitrary cross-section, provided the cross sections 
symmetrical with respect the centroidal axis perpendicular the 
plane flexture. The expressions for f(h2) and change, the core 
radius takes the place the dimension h/6 and cannot expressed 
Eq. (6). But always continuous function and and the 
two limiting cases considered above takes the constant values and 
all along the beam. coincides again with the reduced modulus von 
Karman for the cross-section shape considered. 


Evaluation the Critical Buckling Load 


well known that for less than the only possible solution 
the boundary value problem (9), (10) v(z) while for Py, to- 
gether with the solution v(z) two additional solutions v(z) are 
possible. These additional solutions are symmetrical with respect 
the axis. each value interval having lower limit 
there corresponds only one couple solutions,'+ v(z). Moreover, 
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parameter Hence, the value corresponds the first bifurca- 
tion Eq. (9) and, once the existence this point established, its 
location may obtained follows. 

Let call the solution Eq. (9) corresponding value 
near and higher than Pg. The solution approaches zero ap- 
proaches But since bifurcation point, Eq. (9) also admits 


z=0 
(11) 


z=0 (12) 


Neglecting infinitesimals higher order Eq. (12) may written as: 


where the partial derivatives are evaluated Subtracting Eq. (11) 
from Eq. (13) and neglecting the term obtain: 


tion which, will seen, known. Thus Eq. (14) becomes 
the limit: 


z=0 


Since bifurcation point Eq. (9), Eq. (15) must admit 
equation: 


with the boundary conditions: 


(17) 


The value taken approaches depends the dis- 
tribution the stresses the cross-sections the beam when the 
beam starts buckle. the first case considered above, and 
all along the beam; hence: 


and Eq. (16) becomes: 
(18) 


The first eigen-value Eq. (18) corresponding the boundary condi- 
tions (17) known equal hence: 


531-6 


(19) 


This the tangent modulus load Engesser. 
the second limiting case, instead, and all along the 
beam; hence Eq. (16) becomes: 


(20) 


(21) 


This the reduced modulus load von Karman. 

well known theorem the eigen-values linear allows 
conclude that the nth eigen-value Eq. (9) bound the two nth 
eigen-values the two equations obtained from Eq. (9) substitution 
with and respectively. Hence the buckling load always 
bound below and from above. 

Since the value which buckling start falls between the 
Engesser and the von Karman loads, depending the initial distribution 
the stresses which arbitrary, engineers should base their de- 
signs the conservative assumption These conclusions are 
identical with the conclusions reached Shanley for the case sim- 
plified beam model. 


Graphical Representation the Buckling Phenomenon. 


The graph Fig. gives the deflection v*of the middle section 
the beam function The curve v*(P) symmetrical with re- 
spect the axis and crosses this axis arbitrary point be- 


and hence: 


v*>0 v*>0 


from which 


and 
4 
lyk — ° 
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the radius the section core. instead and the 
slope zero. 


Although the results obtained Mr. Franciosi this paper coin- 
cide with those obtained Mr. Pfliiger the reference listed 
footnote (2), was felt the reviewers the Engineering Mechanics 
Division that the original mathematical approach used Mr. Franciosi 
justifies making his results known the English language the mem- 
bers the Society. 
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